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Abstract 

We calculate the expectation value of the coincident product of two field strength tensors at 



^ \ two loop order in scalar electrodynamics on de Sitter background. The result agrees with the 
J^T ■ stochastic formulation which we have developed in a companion paper [2] for the nonperturbative 
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resummation of leading logarithms of the scale factor. When combined with a previous compu- 
tation of scalar bilinears our current result also gives the two loop stress-energy tensor for 
inflationary scalar electrodynamics. This shows a secular decrease in the vacuum energy which 
derives from the vacuum polarization induced by the inflationary production of charged scalars. 



I. INTRODUCTION 

This is the third and final paper in a series with the goal of establishing a stochastic theory of 
inflation for scalar quantum electrodynamics (SQED). Stochastic SQED is outlined in Ref. 
In this paper and in Ref. [l| we compare the predictions of stochastic SQED with explicit and 
exact perturbative calculations. In particular, Ref. [l| gives the mostly scalar contribution to 
the two-loop stress energy, and here we calculate the mostly photon contribution. 

Stochastic inflationary theory was developed primarily in order to capture the ef- 

fects of small scale stress-energy fluctuations (large momenta) onto large scale fluctuations, which 



are observed by current large scale structure and cosmological microwave background radiation 
measurements. Furthermore, stochastic inflation has also been used to study global structure 
and dynamics [7, Q] of inflationary universe models. Even though corrections to stochastic in- 
flationary theories have not yet been systematically studied, these type of studies are feasible 
within the in-in path integral formalism 9f. A phase space formulation of stochastic inflation- 
ary theory is given in Ref. [1Q[, where a quantum Liouville equation for the Wigner function 
was obtained, which is the suitable quantum generalization of the Fokker-Planck equation for 
the (classical) density function. More recently, a phase space formulation has been used for 
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computing the nongaussianity of cosmological perturbations [111 , 
inflationary models. 

We are interested in the stochastic theory of inflation mainly because it captures correctly the 
leading logarithmic behavior (in the scale factor) of quantum field theories during inflation. In 
Ref. 2j we have proven that a stochastic theory of scalar electrodynamics (obtained by integrating 
out the photon field and by subsequently stochasticizing the resulting effective scalar field theory) 
correctly reproduces the leading secular logarithms to all orders in perturbation theory. The 
agreement is made explicit by a comparison between the mostly scalar contribution to the two 
loop stress-energy tensor in de Sitter inflation with the corresponding stochastic prediction. 
The stochastic theory of inflation goes far beyond perturbation theory .in that it prov.des an 
effective means for resummation of the leading logarithm contributions |2| . Unlike full quantum 
field theories, which cannot be generally solved, stochastic theories are classical, and hence they 
are relatively easily solved. Indeed, solving a stochastic theory is tantamount to solving the 
corresponding Fokker-Planck equation for the density function 5|. 

One of the important application of stochastic theories is that they permit to address the 



question of gravitational back-reaction in accelerating space-times [171 1 1 8f | , and can thus answer 
one of the fundamental unsolved problems of inflationary cosmology: how do quantum effects 
change the dynamics of space-time during inflation. More broadly, this will teach us something 
about the infrared sector of quantum gravity and the cosmological constant problem. Ultimately 
one would like to understand how quantum fluctuations of gravitons affect inflation. The hope 
is that this question can be answered using a stochastic theory of inflation 191 ]. 

Quite a lot is already known about quantum (radiative) effects in inflationary space-times. 
Perturbative analysis of scalar electrodynamics during de Sitter inflation shows that the photon 



acquires a mass during inflation 
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251 ] , which can be of relevance for generation 
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small 
light 



of the cosmological magnetic fields 20|, ]23|, l25|, ]26|, ]27|. The scalar field of SQED also acquires a 
mass during inflation, but — unlike the photon mass 2( — the scalar mass remains perturbatively 



. Similarly, the fermions 28|, [29| of Yukawa theory acquire mass while the scalars remain 



311 ] . In self-interacting scalar theories, the quantum loop effects during inflation can 



contribute to the stress-energy like a phantom field 32], ]33j], whose equation of state is even more 



negative than the pressure of a vacuum filled with a positive cosmological term. Alt houg h this 



generally betokens instability 



34|, the phenomenon in this case is self-limiting |35l . |36| |. The 



quantum effects of g ravitons on the dynamics of fermions and scalars during inflation have also 



been studied 
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39|, 



40] , and it is now known that gravity induces a logarithmic correction to 



the fermionic wave function at one loop while no such correction occurs to a massless, minimally 
coupled scalar. The question of how large quantum corrections to the inflationary observables 
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ly are, has been addressed by a number of authors 
521 ]. but as yet no definite answer can be provided to that question. 
This paper is organized as follows. In section II we introduce the relevant scalar and gauge 
field propagators for de Sitter space and discuss some of their properties. In section III we 
compute the one loop expectation value of the SQED stress-energy tensor in de Sitter space. In 
section IV we discuss some general properties of the stress-energy tensor and its constituents. 
We also employ previous work |l| to give the two loop contributions from the scalar kinetic and 
potential energies. Section V is the central part of the paper, where we compute the two loop 
gauge field strength bilinear. In section VI we obtain the total, two loop stress-energy tensor 
and we discuss the result. Appendices A, B and C give identities and integrals useful to readers 
who wish to reproduce the calculation. 



II. THE PROPAGATORS IN DE SITTER SPACE 



In this paper we approximate inflationary space-times by de Sitter space. The de Sitter metric 
with flat spatial sections is, 

ds 2 = a 2 (r]) (-drf + dx 2 ) , (1) 

where a = a(j]) is the scale factor, which in de Sitter space has the form, 
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Here H is the Hubble parameter and r\ denotes conformal time. Furthermore, it is convenient to 
define a de Sitter invariant function y = y(x;x') between two space-time points x M and x' v . In 
flat coordinates (pQ), 



y = aa'H Ax 



Ax 2 



i\v-v'\ 



ie) 2 + \\x — x'\\ 2 



(3) 



where a = a{rj) and a' = a(r]'). There is a simple relation between y and the geodesic distance 
£ = £(x; x 1 ) in de Sitter space, 



y = 4 sm -— 



(4) 



We work here in the Schwinger-Keldysh formalism |53l . l54j . l55l . |56( |. which is suitable for 
describing quantum field theory dynamics. In short, the vertices of Feynman diagrams have two 
polarities: plus (+) and minus (—). Because pairs of vertices have four possible polarities there 

are four propagators which we denote by iA ++ , iA + _, zA [_ and iA , respectively. They are 

obtained from the Feynman propagator iA = iA ++ by replacing y = y ++ = y(x + ; x' + ) in Eq. ([3]) 
by y+-, y~+ and y , respectively, defined by, 
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A. The massless minimally coupled scalar propagator 



The Feynman propagator for a minimally coupled massless scalar field in de Sitter space obeys 
the differential equation, 

v^DiA(x; x') = iS D (x - x') , (8) 

where □ = (— g)~ l l 2 d^ u yf~^gd v denotes the scalar d'Alembertian, g = det[g Ml/ ], and S D is the 
-D- dimensional Dirac delta function. The addition of a mass term to (18]) results in a de Sitter 



invariant propagator which was first constructed by Chernikov and Tagirov 571 ] ■ However, there 
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is no de Sitter invariant solution in the massless case 

The problem can be circumvented by adding de Sitter breaking contributions to the propaga- 
tor. A natural way of obtaining the form of the de Sitter breaking term is to consider quasi-de 
Sitter space in which the de Sitter symmetry is mildly broken, such that e = —H/H 2 <C 1 and 
e = (here H = dH/dt). One can show 59( that the scalar propagator in such a space can be 



represented as a function of the de Sitter invariant distance suitably rescaled by powers of the 
two scale factors. Remarkably, in the de Sitter limit when e — ► 0, one recovers the massless min- 
imally coupled scalar propagator in D space-time dimensions [33J plus an infinite e-dependent 
constant. Ignoring the constant, the propagator in de Sitter space can be written in the form, 

iA(x;x') = A(y) + kln(aa) (9) 
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Note that the function A = A(y) in Eq. ( TTUl) obeys the equation, 

(4 - j/)yA"(y) + L>(2 - y)A\y) — (D — l)k . (12) 

This implies a relation which will have great importance for this computation, 

K'(v) 

(2-y)A'(y)-k = -^l, (13) 

where we define K(y) as, 

K(y) = (4-y)yA'(y) + (2-y)k. (14) 
The expansion of K(y) around y — is, 



^;4^|ffl'-* + «ML +0((D _ 4m . (15) 



(4vr) D / 2 V 2 / 1 V4/ 2 1 - y/4 

Note that in D = 4 the function K(y) terminates at y 1 ^ / 2 . In fact, up to the terms that are 
suppressed as (D — 4), K(y) is just —2{D — 2) times the conformal scalar propagator. Another 



5 



important relation is, 



(4 - y)yK"{y) + D{2 - y)K'{y) = (D - 2)K(y) 



(16) 



Also note that dimensional regularization implies the following coincidence limits (y 

H D ~ 2 /tt£>\ , / D 
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B. The vector propagator 



The Lorentz gauge vector propagator in de Sitter space was discussed by Allen and Jacob- 



son 



60l | . and a minor error in their analysis was recently corrected 6JJ. The propagator obeys 



the equation 



g (Ug pl * - R pp ) i ^A u (x; x') = 8 p iS D (x - x') + ^d p d'jA(x; x' 



(20) 



where □ denotes the vector d'Alembertian, R pp is the Ricci tensor and iA(x;x') represents the 
scalar propagator ([9]). The transversality of Eq. (1201) follows from the Lorentz gauge condition 



Ay 



(x; x') = d p g plx v 73 ^ i (x; x 1 ) = 
(x]x) = d! <jg'° v \[-g'i Au (x;x') = 



(21) 



where g ni = g pp (x), g' uv = g cru (x'), and V M is the covariant derivative operator. 

On de Sitter background the photon propagator can be written in the following de Sitter 
invariant form 



{x-x') = B{ y y 



d 2 y 



C{y) 



dy dy 



(22) 



dx^dx' v dxt 1 dx w 

We can use the Lorentz gauge condition (j2ip to express B(y) and C(y) in terms of a single 
function j(y) 
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A(D-1)H 2 
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4(D-l)H 2 
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{2-y)i{y)-{D-l) 1 {y)\ . 



(23) 
(24) 



The propagator equation will follow provided 7(1/) has the right singularity at y = and obeys 
the differential equation, 



(4-y)y>y"(y) + (D+2)(2-y)j'(y) - 2(D-l)j(y) = (D-l) (2-y)A\y) - k 



The unique solution for j(y) is [61| 
D-l H°- 2 
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An important simplification for the vector propagator is "the outer leg identity" 

1 dy d 2 y \ fn \ At , \ , ] K'(y) dy d 2 y 



[x: x 



AH 2 dx^ dx^dx'p 



(2-y)A'(y)-k 



(25) 



A(D-2)H 2 dx^dx^dx 



(26) 



(27) 



Our convention throughout this paper is that square bracketed indices are anti-symmetrized: 
= (H llu —H UfJ )/2. The outer leg identity can be easily obtained by making use of Eqs. (I22T 
251) . Note that Eq. ( 1271) allows us to express the integrals entirely in terms of the A-type 
propagator (IU HTT1) . without involving the complicated function 7(1/). This is a very welcome 
simplification because the series for A(y) terminates in D = 4, whereas that for 7(7/) does not. 



III. THE ONE LOOP STRESS-ENERGY TENSOR FOR SQED 
A. The SQED action 

When expressed in terms of the renormalized fields, the action of scalar electrodynamics is, 

SsQED = J d D xCsQED (28) 

Csqed = - (1 + 6Z 2 ) {V^Y (V v( p) g^V=9 ~ \{± + SZ^F^F^g™^ (29) 

Here V>^ = + ieA^ is the covariant derivative operator, A^ is the photon field, <p is the scalar 
field, g = det^jy] is the determinant of the metric tensor g^ u , and g pa is its inverse, g pp g pu = 8^ v . 
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Note that we have chosen the renormalized values of the conformal coupling, the scalar mass 
and the scalar self-coupling to vanish. There is no need for a mass counterterm because mass is 
multiplicatively renormalized in dimensional regularization. However, full renormalization does 
require a conformal counterterm and a scalar self-coupling counterterm. The various one loop 
counterterms are 12J, 122J , 
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-V»p-l)+V>(2)+v(y-l)-V(2 




(33) 



Note that, while the finite parts of the two field strength renormalizations are arbitrary, the finite 
parts of 5£ and SX are fixed (at one loop order) by requiring that the scalar effective potential 



be as flat as possible 



2]- 



From the action ([281 one obtains the stress-energy tensor in the standard manner, T a p = 
— (2/ \ / —g)8SsQED/fig a ^- Its expectation value in the presence of Bunch-Davies vacuum \ fl) can 
be broken up into a part derived from the Maxwell field strength and a part derived from the 
scalar kinetic and potential energies, 



(n\T aP (x)\n) = (fi\T al3 (x)\n) M ^ wcll + (n\T aP {x)\n) 



scalar 



Maxwell 

(fi|T a/3 (x)|f2) sca i ar 



SZ^^F^F^x)^) 



<W - -gapg^j 2(1 + 5Z 2 )(n\{V„ip{x))\V v ip{x))\V) 
+ 2<5£ (G af3 + gapg^VpV,, - V Q V^) (Q\^(x)^(x)\n) 



-SXg a p(QMx)^(x)] 2 \Q) . 



(34) 



Here G a/ 3 = R a p — \g a pR is the Einstein curvature tensor, and denotes the covariant deriva- 
tive. On de Sitter background the various curvatures become, 



G af3 = ~(D-l)(D-2)H 2 g aP 



R a(3 = (D-l)H 2 g aP and R = D(D-1)H 2 . (35) 
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B. One loop stress-energy tensor 

At one loop order the various counterterms are irrelevant and we have, 

(Q\T a(3 \Q) lloop = [g^ggg va -\g a ^ p g va ){^\F^F pa {x)p) lloov 

+ (s 8 V-^) 2<n| W fl ^- (36) 

The electromagnetic contribution derives from the coincident product of two field strengths, 
(n\F^(x)F p<7 (x)\tt) lloop = (dM„A a ](x;x')) + fe^z[ M A p ](x; x') 

\ J x=x' V 

- (d„d' p ilA a ](x;x')) -(d lt d' a i[ v A p ](x;x , j) , (37) 

V / x=x' V / x=x' 

while the scalar contribution comes from the scalar kinetic energy, 

(n\d^*(x)dMx)\tyn oop = (d p d' u tA(x;x')) . (38) 

V / x=x' 

Making use of the outer leg identity ff27|) we find that the nonvanishing terms are 
(QlF^F^xm, loop = ^— {(fffay) (dy&y) K'iv)}^ = g p[ ^ ]a . 



(39) 



When this is inserted into Eq. (jT9~I) one arrives at, 



jjd y(Fj 1) 

(n\F pu (x)F pa (x)\n) lloop = - — -TT—r- — ^(g^gua - g^g vp ) , (40) 

(47r)T F ^ + 1 j 

Similarly, we make use of relation ([TBI for the scalar propagator (Q to find, 

H D T(D) 

i loop _d / \ g^y ■ 

(41) 

(47r)T 2 r(f + lj 

The one- loop stress-energy tensor is easily obtained from Eqs. ( l36l) . ( 1401) and ( 14T1) . The 
electromagnetic contribution is 

(n\r m\ - HDv{D) D ~ A n ^ U2) 

\«|-ta/3|i'/l loop, Maxwell — ~ 7 \ : #0/3 ► U, 1 42 J 

(4vr)fr(f + l) 4 

while the scalar contribution reads, 

/OIT IO\ - H ° r (^) " 3if4 ^ 

\* ' | a/3 * '/ 1 loop, scalar — " ~d — J~ n ~ g a /3 > : „ 2 ffa/3 • (4.5 J 

(47r)- r (f + lJ 2 167r 



Summing the two contributions (|42H43|) results in the final expression (1361) . 

(47r)-2rffj lb7r 

These one loop contributions to the stress-energy tensor in de Sitter space are of the form a 
constant times g a p, hence they can be absorbed in the cosmo logical term. 



IV. GENERAL PROPERTIES 



The one loop stress-tensor and its constituents are de Sitter invariant. Their tensor structures 
are comprised entirely from products of the metric, and the coefficients of these products are 
constants. That is not true of higher loop contributions. These receive de Sitter breaking 
contributions from the undifferentiated scalar propagator (Q and from vertex integrations that 
reach back to the finite initial time. Spatial homogeneity and isotropy is preserved but tensor 
structures can involve the unit timelike vector a£°, and coefficients can depend upon ln(a). Of 
course it is the temporal growth implied by factors of In (a) that gives this computation its 
physical interest, and it was to reproduce and sum up the leading powers of these logarithms to 
all orders that the stochastic formulation was developed. 

The expectation value of T^ u at a general order derives from the expectation values of four 
composite operators: -F^F pcr , (D^ipYD^, ip*ip and (ip*<p) 2 . It is useful to parameterize these 
expectation values as follows: 

(QIF^F^Q) = Eg^g a]u + Fa 2 5lg u][ J° p] , (45) 
{Q,\{D^)*D^\Q) ee Jg, v + LaHpl , (46) 

(J%V|n> = ^ , (47) 
5X(Q\{(f*if) 2 \Q) ee N . (48) 

The expectation value of the stress-tensor can be given a similar form, 

(n\T, u \n) = pg, u + (p+p)a 2 5X • (49) 

The coefficient functions E, F, J, L, M, N, p and (p + p) all take the form of H D / (4?r)^ times 
functions of the scale factor and the dimensionless loop counting parameter e 2 H D ~ A / (An)^ . In 
SQED there can be at most one additional factor of ln(a) for every additional loop [2j, and it 
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turns out that the temporal coefficient functions F, L and (p + p) are always down from this 
limit by at least one overall power of ln(a). Note that we have already encountered this pattern 
at one loop where there are no powers of e 2 H D_4 /(47r)^, the coefficient functions E, J and p 
are constant, and the three temporal coefficient functions accordingly vanish. 

In addition to factors of In (a), the various coefficient functions generally depend in a compli- 
cated way upon inverse powers of a that redshift to zero at late times. This can be illustrated by 
the two loop result for the expectation value of the stress-energy tensor of a massless, minimally 



coupled scalar with a quartic self-interaction 32|, |33j. We report this with a minor change of the 



renormalization scheme to bring it into conformity with our convention of not employing a mass 
counterterm, 

5A hn XH 4 J 1 ,. . /15 7 E \, ,, 37 lE 1 ^ n 2 -4 



l P J V = ~ 87tG + (2^{- 8 ln (a) " b + 32 j ln(a) + 144 + 48 " 24^ (n+l^a 
+ 0(\ 2 ), 

( \ ^H 4 f 1 , / a 37 lE 1 1 1 ^n + 2 1 1 ^. x2 . 

l P + P V = (27Fl-^ ln(G) -2^-% + ^^-2l^^I^f + ° (A) - 



n+l 



The larger number of infrared logarithms — two for each extra loop in this model, as opposed 
to one in SQED — is simply explained by having four, as opposed to two, undifferentiated 
scalars in the interaction [2J. What concerns us at the moment is the terms which fall off like 
powers oil /a. These terms derive from the lower limits of conformal time integrations. Because 
they are separately conserved it has been conjectured that these terms can be absorbed into 



perturbative corrections of the initial state 



If so, they are an artifact of the initial state 



rather than a universal feature of inflationary particle production. Because these sorts of terms 
fall off at late times they play no role in the stochastic formalism we are checking. Hence there 
is no point in struggling to retain these terms and we shall not bother to do so. This allows 
vast simplifications from partially integrating without concern for surface terms, which we have 
already done in evaluating the scalar coefficient functions J, L and M at two loop order [l|, 
and which we shall be using throughout this work e.g. when extracting d'Alembertians from the 
inner loop. 

An all-orders relation between the scalar coefficient functions is implied by the operator equa- 
tion Q, 

□faV) = 2g^{D^YD uV + -^L- itVV + t^tt (V?V) 2 • (52) 

l+dZ 2 l+oZ 2 



11 



Taking expectation values and substituting equations (I46H47I) gives, 
- M" - (D-l)M' = 2D J -2L + 2D [ D Zl^ M + 



N 



(53) 



1+5 Z 2 1+5 Z 2 

where a prime denotes differentiation with respect to ln(a). A slight rearrangement produces a 
key result, 



2D 



[l+5Z 2 )J + 5i(D-l)M 



'l+Z, 



M" + (D-1)M' 



2(1 + 5Z 2 )L + N . 



(54) 



It is now time to work out the Maxwell and scalar contributions to p and (p + p) in terms 
of the various coefficient functions. We find the Maxwell contributions from relations ( 1341) and 

ED, 
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:i+5z 3 



p+p) 



Maxwell 



Maxwell 



-\(D-l)(D-A)E + ~(D-3)F 



1 



[1+5Z 3 )-(D-2)F. 



(55) 
(56) 



The scalar contribution to (p + p) is almost as simple, 



p + p) 



scalar 



[l+5Z 2 )2L + 5£ 



-2M" + 2M' 



(57) 



However, it is best to modify the scalar contribution to the pressure using Eqn. fl54|) . 



pj = -{D-2) {1+5Z 2 ) J+8£{D-1)M + (l+5Z 2 )L-25t M"+(D-2)M' - -N , (58) 

scalar L J L J 4 



[l+5Z 2 



D-2 
2D 



M"+{D-l)M' 



-25£ 



M" + (D-2)M' 



.1(1+SZ 2 )L+(^)N .(59) 



This final form fl59|) is quite significant because it precludes ln(a) contributions from the scalar 
at two loop order. This is a key prediction of the stochastic formalism [2J. 
Stress-energy conservation provides an important accuracy check, 



V"{n\TM = -Ha6° u ( p + p )'- v '+(D-l){p+p) 







(60) 



Recall that a prime denotes differentiation with respect to ln(a), and note that the powers of 
l/o we ignore should be separately conserved. The Maxwell and scalar contributions are not 
separately conserved but the Heisenberg equations of motion for SQED relate their divergences 
to the expectation value of a field strength contracted into a current, 



v^niT^in) 



Maxwell 



-V^n|T^|fi) scalar = -ie(l + 5Z 2 )g pa (n F w <p*D a (p—(D a (p)*(p 



n 



(61) 
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Although we shall not give the derivation it is simple to evaluate fIBTl) at the order e 2 relevant to 
this analysis, 



-ie{l + 5Z 2 )g pu (vL 



MP 



-HaS) 



tp*D a tp-(D a ip)*(p 

2o-2£>-4 

-T(D-l) 



n 



(4tt) 



D 



£>-4 



- 10 + 0(a- 1 ,D-A) +0{e 



(62) 



This is probably the right point to discuss finiteness. Conventional renormalization only serves 
to remove divergences from noncoincident 1PI (one particle irreducible) functions. Composite 
operators such as F /ii/ (a;)Fp (7 (x) and the others require additional, composite operator renormal- 
ization to remove their divergences. This is evident from the divergence of expression ( 1621) . One 



renormalizes composite operators by adding a series o 
vergences order-by-order in the loop expansion 



62 



counter-operators which remove the di- 



631 ]. As with conventional renormalization, 



there are ambiguities in finite parts, which are resolved by imposing additional renormalization 
conditions. There is absolutely no reason to burden an already formidable analysis by bothering 
with any of this. The stochastic formalism which we are trying to check provides unambiguous 
predictions for the dimensionally regulated expectation values of composite operators 3], so we 
may as well leave the regulator on and verify these predictions. 

A related point concerns the stress-energy tensor. Unlike the other composite operators, 
its expectation value is a 1PI function of an enlarged theory; specifically, its represents matter 
contributions to the 1-graviton 1PI function in gravity + SQED. It must therefore be that 
divergences in the expectation value of the stress-energy tensor can be absorbed with the addition 
of purely gravitational counterterms which all degenerate on de Sitter background to changes 
in the bare cosmological constant. This was shown at two loop order in Xip 4 theory by explicit 



computation 



32 



331 ]. and the ability to change the result by an additional, finite renormalization 



of the cosmological constant is evident in expression ([50]) . For SQED this means three things: 

• Divergences in the temporal coefficient functions F and L must cancel when they are 
formed to give the full (p + p); 

• The parts of p which contain factors of ln(a) must be finite; and 

• The parts of p which are constant can diverge. 

The first two facts provide more accuracy checks; the third means that there is no point in 
working very hard to determine the constant contributions to p because any constant part of p 
can be subsumed into a renormalization of the cosmological constant. 
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It remains to summarize what is already known about the coefficient functions. The Maxwell 
coefficient functions (|45p take the form, 

H D f2r(£>-l) e 2 H D ~ A 



E 



(4tt)t\ T(f + 1) (4tt)t L 



E 1 \n(a) + E 2 + 0(a- L ) + 0(e 4 ) 



H 



D 



2 TjD-i 



(47r)^ 



+ 



e 2 # 



(47T)" 



-In 



(63) 



(64) 



The point of this paper is to compute the D-dependent numbers E±, E 2 and F±\ They give the 
order e 2 H D ~ A / '(4tt)% electromagnetic contributions to p and (p+p), 



V 



H 



D 



Maxwell 



(47T 



{D-A)T(D) i e 2 F 
"747 



2 rrD-4 



4r(f +i; 



--(D-l)(D-4)^iln(a) 



-~{D-1){D-4)E 2 + i(£>-3)F! + o(a~\ D-A 



2T(D-1) 
Z)(£>-2)(£>-3) 

+ ^(e 4 ) 1,(65) 



2 rrD-4 



+ 



e 2 # 



i(D-2)F 1 + 0(a- 1 ; 



0(e 4 



Maxwell (47r) "2" I (471")" 

Note that Eqn. (162]) implies a relation between the coefficients E\ and F 1; 

i(2?-l)(D-4)^ + \{D-1){D-2)F X = - ^"^ - 20 + Op-4) 



The stochastic formalism predicts [2J 

8r 2 (D-i) 

r( 2 

16 + C( J D-4) . 



Ex 



16r(D-l 



^p-l)+V(2) 



L>-4 

So we can use Eqn. (1671) to infer what Fi should be, 

ierp-1) 



+ 0(L>-4) . 



(66) 

(67) 

(68) 
(69) 

(70) 



(D-l)(D-4) 

No comparable check involves E 2 , and we see from (|65|) that F 2 has no physical relevance because 
it can be absorbed into a renormalization of the cosmological constant. 

The coefficient functions (1461) associated with the scalar kinetic energy are [lj], 



J 



H 



D 



(An) 



hfl(x ^^p- 



2T(f + 1 



r( 



/-D7T 

-2 ln(a) + 7r cot (^ _ ^ _ 



+ 5Z 2 x 



rp) 



2T(f + 1) 



e 2 H D - 4 
(47r)f 



i6r(D-i) 



p-l)(Z>-4) 



17+27r 2 + C(a~ 1 ,L)-4 



+ 0(e 4 ) 



(47T) 



+ 



IjjD- 



(4tt; 



ior(D-i; 



(D-l)(D-4) 3 



- + 0(a _1 ,D-4 



+ C(e 4 ) 



(71) 
(72) 
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Recall that 5Z 2 and <5£ are of order e 2 H D 4 / (An) 2 and are given in expressions (I5U1) and 
respectively. The coefficient function M was defined in Eqn. (1471) and takes the form [1], 



M 



H D \Y(D-l) 



(4tt)t 



r 



D N 



/7T.D 

2 ln(a) — 7r cot (^ _ ^ _ 



<5Zo x 



r(f) 



-21n(a) + 7rcot(^— J 



+ 



e 2 H 



(An) 



-108r(D-l) 
( J D-l) 2 (D-4) 



-40- 



16vr z 



ln(a) +Mi + C(a" 1 ,L'-4 



+ 0(e*) . (73) 



The constant Mi has not been determined but it drops out because only derivatives of M enter 
expression (159!) for the stress-energy tensor. Note that even the lowest order contribution to M 
contains an infrared logarithm, so the t loop result can have up to I factors of ln(a), as opposed 
to only I — 1 in J. This is entirely in conformity with the rule that each additional loop brings 
at most one additional factor of ln(a) [2|. The final coefficient function ( T48l) is, 



N 



H D \ , 5XH D ~ i 2T 2 {D-l) 



(4tt)t 



(4tt 



r 2 (f) 



2 ln(a) — 71 cot 



0(e e 



(74) 



Because SX in (I33p is of order [e 2 H D ~ A /(Atc)%] 2 , the function N does not affect the two loop 
stress-energy tensor. 

We infer the scalar contribution to the stress-energy tensor by substituting these coefficient 
functions into expressions (1591) and (1571) . 



scalar (ATT) 2 



H D J (D-2)T(D) e 2 H 



2 rrD-4 



2T(f + 1) 



(4tt 



x ln(a) + 



94 

y 



2r(z>-i) 

(D-l)(£>-4) 
An 2 + o(a- 1 ) D-A 



P + P 



scalar 



H e 2 H D ~ 4 
+ 

(4tt)t 



(4tt) 



8r(D-l) 20 
(D-l)(D-4) + T 



+ C(a~ 1 , J D-4 



+ ^(e 4 )^, (75) 
0(e 4 U. (76) 



Note that there are no infrared logarithms in (l75j) at this order. Note also that ( |76j) obeys the 
partial conservation identity (l61H62p . which represents an additional check of the two loop scalar 
stress-energy calculation of Ref . H. 



V. THE TWO LOOP GAUGE FIELD STRENGTH BILINEAR 



Three diagrams contribute to the VEV at two loop order. They are depicted in Fig. 1. 
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FIG. 1: Two loop contributions to (^l\F fll/ (x)F pcr (x)\n) . 
A. The counterterm diagram 

The diagram on the far right of Fig. 1 derives from the photon field strength renormalization. 
Covariance and index symmetries dictate its form, 



l r . t . = 4i8Zj d D x'{d {lx i 



(x;a; / )v / -F 



g'a/3 fj5_ faS ifij 



,x: x 



(77) 



where stands for the covariant derivative operator with respect to x >tl . To evaluate this 
diagram we first interchange the order of covariant differentiation in the term proportional to 
gia&gip^i an( j tn en use the Lorentz gauge condition to get, 



a}Ap 



[x\ x 

+± 



<r] A/3 



x; x 
-i +± 



Now we use the appropriate Schwinger-Keldysh generalizations of the propagator Eq. (T2UD . 

i5 a J D (x-x') + v^g' al3 d' p d a iA ++ (x-x') , 
(x; x') = . 



V-9' 




i 


*Ap 


V-9' 




% 


ct Ag 



(78) 

(79) 
(80) 



The anti-symmetrized derivative with respect to x p makes the second term in (179]) drop out, so 
the result for the far right diagram (J77J of Fig. 1 is (cf. Eqs. f l37j) . fl39l) and fHDl). 



l c . t =-45Z 3 Jd D x'd [p i 
= —5Z 3 lim ld,,d'i 



i/iAq, 



3 imi w/n/.. 



25Z-.H 2 



2A'(0) - k 



(z; x')d [p 5 a a] 5 D {x-x') 

++ 

(x\x')-d v & p i ^A^ {x]x')+d v d' a i 

9pp9va 9pu9up 



(x]x')-d p d' a i U A P (x;x') 



-5Z 3 x(n\F^(x)F pa (x)\n) 



one loop • 



From expression fl55|) we see that the contribution of X c . t . to the stress-energy tensor is exactly 
canceled by the contribution from SZ 3 times the one loop result. 
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B. The figure 8 diagram 



The middle diagram of Fig. 1 is, 



Zoo = -8*e 2 / d D x'y^\ d^i L]A J Or; x')d [p i L A J (x; x') - (+-) \g' al3 iA(x'; x') . (82) 



Making use of the outer leg identity (12TI) transforms this to, 



T 



— 



2(L>-2) 2 if 4 



A'^A^;!') -(+-) 



<9y <9 2 y 



to/3 



<9x dx 1 ^ <9x' Q <9x^ <9x [ CT <9x^ 



Now use the contraction identity 



65j, 



^ j a0 &*y 



dx u dx' a ^ dx >l3 dx a dx v dx a 

The final term drops out owing to the anti-symmetrizations, so the middle diagram becomes 



;4) 



T 



f d D x , ^/^iA{x , ;x')^ ; g v][(T d!l 



(D-2) 2 

-2ze 2 
(D-2) 



dx^ ulla dxrt 




K'(y ++ ) 



1 2 



^£(±) I d D x'^iA(x' ] x') x ^K{y +± )g v]W ^K{y +±J 



Q x p\ 



55) 



From Eqs. (TroTTTTT) and (flT|) we can read off the coincidence limit of the scalar propagator, 



iA(x';x') = A(0) + 2k ln(a') 



#/j-2 r(rj-i) 
(47r)f r(f) 



C0S (f i, ) r lij r l 1 -Yl + 21n ( a '' 



(86) 

We will see that the leading logarithm contribution E\ derives entirely from the 2fcln(a') term 
in the integral (I8~5l) . 



C. The two leg diagram 

The far left diagram on Fig. 1 is the most difficult in any gauge. To simplify it we must derive 
new identities for contracting adjacent basis tensors into one another. Let us first make the 
following definitions for the three coordinate separations and their associated length functions, 

Ax M = x M -x' M y = y(x; x') = aa'H 2 Ax 2 , (87) 
Ax'» = x'^-x"" y' = y(x'; x") = a'a"H 2 Ax' 2 , (88) 
Ax"» = x^-x" M =^ y" = y(x; x") = aa"H 2 Ax" 2 . (89) 
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The three fundamental "adjacent" contractions we require are, 



dV 9' P °^L = H 2 {2y + 2y'-yy'-2 y "} 



dx'p dx' a 
d 2 y ,„ dy' 



dx^dx'P 9 dx' a 
9 2 V ,„ d 2 y' 



H 2 {(2-y 



^ dy 2 dy" \ 
dx^ dxP J 



:9 



-2H 



2 d 2 y" 



H 



2 dy dy' 



dx^dx'P dx' a dx" v dx^dx" v dxP dx" u 

For x"p — > x^ these relations degenerate to the well known identities 

dy ,pc dy 



dx' p 9 dx" 7 
d 2 y _, pa dy 



(A-y)yH 2 



;9 



dxP-dx'p dx' a 



H 2 (2-y) 



dy 
dxP 



d 2 y .op d 2 y 



(90) 
(91) 
(92) 

(93) 
(94) 
(95) 



dx v dx' a dx'Pdx a dx v dx a 

Note that the last of these is just the contraction identity (1541 . The contractions on the "other" 
site x" are trivial to obtain from these by interchanging x lfi and which interchanges y and 
y", while leaving y' unchanged. Two other useful identities are, 

„ d 



: y = H I {2-y)g i 



and 



Uy = D{2-y)H 2 



v dx 11 ' 

where V u is the covariant derivative and □ is the scalar d'Alembertian. 
We are now ready to tackle the leftmost (two leg) diagram of Fig. 1, 

J = 8e 2 ^(±)(±) J d D x'^'d^i[ v] A^{x-,x') j d D x"y^d [p t 



(96) 



CT ]A 7 



[X] x 



X 



g >*P g »is) d'iA ±± (x'; x")%iA ±± (x'; x") - iA ±± (x'; x")d'd'^A ±± (x'; x" 



(97) 



The Lorentz gauge condition (12T]) allows us to partially integrate the dp so that the inner loop 
consists of just a single term, 



J = 16e 2 ]T(±)(±) / d D x'^8\ 



u]A a 



(x;x') / d D x"^g"d [p i u] Ary (x;x") 

J +± 



xg' aP g"^dUA ±± (x'; x")d%iA ±± (x'; x") + 0{aT l ) . (9£ 



The surface terms we have neglected can only involve the initial time surface in the Schwinger- 
Keldysh formalism 64|. As explained in section III, these surface terms fall off like powers of 1/a 
and play no role in the check we are making of the stochastic formalism. We shall consistently 
ignore them. 
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Expression (1981) is already a significant simplification because it means we do not have to worry 
about either the undifferentiated logarithm term or the local contribution. It also simplifies the 
tensor algebra. The next step is to act with the derivatives on the two inner loop propagators 
and expand, 

^A ±± (*V')c^A ±± (*V') = A' 2 (y')^^ + kHa'5U'(y^ ' V 



dx'Pdx" 8 p Kyj dx" 5 

+ kHa"5° s A'(y')-^ + k 2 H 2 a'a"5p° 5 . (99) 

The term proportional to k 2 drops out in the polarity sum, so we really only have three terms 
from the inner loop. 

Let us define I[f(y)] as the indefinite integral of f(y) with respect to y. By making use of the 
identity, 

9 ' 9 , 1 2 [A'(y') 2 } = ^L R ^- K A\y') 2 + ^ -J [A'(y'f] , (100) 
dx'P dx" s L yy ' J dx' 13 dx" s Ky ' dx'^dx" 5 K J 

we can rewrite Eq. (1991) as 

d'^A ±± (x';x")d^A ±± (x';x") = -^- I[A ^] + ^^W)] 

+ ^ + ^ [warns AW] 

+ ^^^ 2 Ma')lnK)] • (101) 

Observe that when Eq. (110 II) is inserted into Eq. fl98|) . all terms but the first vanish upon partial 
integration as a consequence of the Lorentz condition (EH), again up to 1/a surface terms. With 
this Eq. (19"5|) becomes 

2 = ^pi E^W / ^'V=?^W) / d D x"^K'(y'l ± ) 



x 



^ d "y a >aP gV Q ^7 ^! gV x jr4«f„' 11 + Ofa- 1 ) (102) 
dx^ dx^dx'" 9 dx'^dx"* 9 dx^dx^dx"^ [ W±±)l^^ a )-^) 



where we used the outer leg identity (1271) twice. 

Note that one can add an integration constant to the first term in (11011) : I[A' 2 (y')] — > 
I[A' 2 (y')} + c/, provided the appropriate change is also made in the second term, I 2 [A' 2 (y')] — > 
I 2 [A' 2 (y')] + ciy'. We have just seen that I 2 [A' 2 (y')] in ( 110 II) can be partially integrated to zero 
as a consequence of the Lorentz condition (121]) . up to 1/a surface terms. Hence it must be that 
any constant term in I[A' 2 (y')] must fall off like powers of 1/a. 
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Equation (11021) can be further simplified by making use of the identities (1911) and (184")) and 
the antisymmetry in fi and v and in p and a to obtain, 

,2 zj2 



x ^^ +± )^[^^:j^^ 2 (^ ± )]+0(^ 1 )- (103) 

We shall often encounter integrals of this general form, 

£(±)(±) / ^y/=P^K(y) j d D x"^f^K{y") x F{y>) . (104) 

Our basic strategy for evaluating them is to extract d'Alembertians for the y' dependent terms 
using the procedure described in Appendix A. In addition to d'Alembertians, this procedure will 
generally also produce a delta function and a finite remainder for which D can be set to four, 

F(y') = \G(y')] + const, x 5 D (x' - x") + R{y') . (105) 



H 2 . 

We partially integrate the □' to act upon K(y), ignoring the order 1/a surface terms. To see 
the effect of acting □' on K(y), note first that, up to delta function contributions arising from 
factors of y 1_D//2 , 

S/(y) = (4 - y)yf\y) + D(2 - y)f(y) . (106) 
In view of relation ( TLB"]) , and the coefficient of y 1 "^ in the expansion i fTBT) of K(y), we have, 

Jf2 K (y ++ ) = -^^-iS D (x-x') + (D-2)K(y) and ^K(y + _) = (D-2)K(y) . (107) 

Hence the general reduction ( 11051) will produce four sorts of terms: 

1. A less divergent, 2- vertex integral from the factor of (D — 2)K(y) in fjlOTjl ; 

2. A potentially divergent, 1-vertex integral from the delta function in fll07[) : 

3. A potentially divergent, 1-vertex integral from the delta function in (j!05P ; and 

4. A finite, 2- vertex integral from the remainder R(y') in fl 105ft . 

Our strategy is to extract another d'Alembertian from the type-1 terms and continue in this 
way until the limit D = 4 can be taken in all 2-vertex integrals. Appendix C explains how to 
evaluate these. It turns out that the 1-vertex integrals of type-3 can be usefully combined with 
the Figure 8 diagram (185]) to cancel an overlapping divergence. 
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Note that the integral (11031) is symmetric under the exchange of the inner legs, x' < 
which y <->• y" and y' <->• y' '. Because the integral is quartically divergent, it suffices 
three most divergent terms in the expansion of I[A' 2 (y' ±± )}. From 



A'(y) = 
we find 



H 



D-2 



1 ( ~ I> + f + l) ( y 

(47r) D / 2 4| ^ T(n + 2) 



-f+i °° 



v ^ r (n + D - 1) /y 
S r(n + f) U 



n— 1 



-> x", under 
to keep the 

(108) 



H 



2D-4 



r 2 (f) 



(4vr) s 



2/' 



1-D 



+ 



2/' 



D-l \ 4/ L> - 2 V 4 

3-D 



2-D 



+ 



+ 



D(D + 1) fy> 



AY{D) 



2(D - 3) V4 



D(D + l)(D + 2) fy' 
6(£>-4) 



r(f)r(f+i) (d-2) V4 

4 ~ D d+6 r(D + i) 
"^4r(f)r(f +2 ) 



x 2 



" 2 



, /-o ( /5-D /3-D/2 /A I 

+ o[y ,y , y ) > ■ 



(109) 



Note that the terms of the order y'^ u ' , y' z Uj ' L or higher in Eq. (11091) are suppressed by (D — 4) 1 . 



/ 4 " D , y' 2 ~ D/2 

We are now ready to extract d'Alembertians from I[A'(y') 2 ] defined in Eq. (11091) . Making use of 
relations in Appendix A we obtain, 

2-D 



I[A'{yT 



H 1D-A V 2 (|) f U , 

' (4tt) d 4 | H 2 



+ 



y_ 

(D — l)(D — 2) 2 V4 
2p + l) 



+ 



(D-l)(D-3)(D-4) I V 4 
2(£> + 1) ( D I 



3-D 



(D-4)(D-6) \2(D-3) (D-l)(D-4) 
2 { D(D-2)(D + 1) 4 



l-4r 



4-D 



(110) 



r(i>) 



+ 



D - 4 ^2(D - 1)(£> - 3)(D - 4) D - 2 T (f ) T (f + l) 

(D - 1)(D - 3)(Z> - 4) T (§ - 1) {Ha>)D 1 J 
~£>(£> + l)(L> + 2) 6(D + 1)/ 4 D 



+ 



6(D-4) £>-6 

l>+6 r(D+i) 
^4r(f)r(f+2) 

L> + 2 /" L>(L>+ !)(£>- 2) 



+ 



(D-l)(£)-4) 2(D-3) )\ V4 



4-D 



2-~ 



+ 



T(D) 



2 I 2(D-l)(D-3)(D-4) D-2T (f ) T (f +l) 
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When inserted into the integral fll03j) the last three lines in Eq. f II 1 j) result in the finite, 
type-4 integral, 



Zfin.l — — e 



2 H 2D - 2 JD 
2(Aix) D 





±± 



l)]T(±)(±) / d D x'a' D d D x"a" D 







x 



101n(^) +0(a~\D-A) 



(111) 



where we dropped the order D — 4 terms as well as an irrelevant constant. The delta function 
in (11 101) leads to the type-3 integral, 

^■ = - fe2 (^(^fe^SS ? (±) / ^' D ^^*^K^) . (H2) 



The d'Alembertian terms Eq. ( 11 101) produce a type-1 integral, 

jy2L>-2 



2k 



(4vr) D 



r (!) r (I" 1 ) g (±)(±) / / ^"v^ 



+ 



(£> — 1){D — 2) 2 V4 



2-D 



2(D + 1) 



(D — 1){D — 3)(D — 4) I V 4 



2(D+1) 



+ 



3-D 



D 



(D-A)(D-Q) [ (D-1)(D-A) 2(D-3) 



4-£) 



(113) 



P-4) 



£>(£>-2)(D+l) 



2(D-l)p-3)(£>-4) D-2r(f)T(f + 1) 



2-J 



and a type-2 integral, 



,H 



2D-A 



sg,2 



(47T) 



D 



2T 



D 



D 



d 



r ( T- 1 ) E(±) / ^V = ^*(v + ±)M* 



X 



5 



4 



(Z>- l)(L>-2) 2 V4/ 
2(D+1) 



2p + l) 



(D - 1)(D - 3)(D -4) V^4/ 



3-D 



4 



(D-4)( J D-6) 



+ 



(D-l)(D-4) 2(D-3) 



4-D 



(114) 



D-4 



D(D-2)(D+1) 



r(z>) 



2(D-l)(D-3)(D-4) D-2r(f)r(f + 1) 



2/V - 
4. 



We have seen that the total integral f 1 1 3 j) breaks up to the sum, 

Z = 2fln,l + 2gg,l + ^sg,2 + Xi + C(a x ) . 



(115) 
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The next step is to extract another d'Alembertian from the y' dependent terms of Eq. (11131) . 
The result is: 



(4vr; 



x —K(y +± )g v][a —^K(y , l ± ){ ^ 



Q x \p, 



±± 



+ 



-5' 



{D-l){D-2) 2 {D-S){D-4) \ V 4 
4 



3-D 



4 D + l\ y' 



(£>-l)(.D-4) 2 (.D-6) V P-2) 2 D-3; V4 

t\ 2 



4 

4-£> 



_D - 
_> 



(D-l)(£>-3)(£>-4) 2 2 ^ + 1 ^) (4 



(47r) D / 2 



(D-1)(D-2)(D-3)(D-4)V(§) (a'tf) 



i5 D (x' - x") 
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+ 



86 



3( D -4)-9( D -4)-|- 51n (l)^ ln2 (ll + °( B - 4 ) 



■ (116) 



where we expanded the terms in the last line in powers of D — 4. Just like I in Eq. (11151) . I\ 
can be decomposed into four sorts of integrals, 
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The contributions of type-4, type-3 and type-2 are, respectively, 
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and 
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The type-1 integral is, 
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This time there is no type-3 integral, 

1l = 2sg,5 + Xfi n ,3 + 2fin,4 + 0(d *) . 



(122) 



The contributions of type-4, type-1 (with D = 4 taken) and type-3 are, respectively, 
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and 
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The finite, 2- vertex integrals X firijl , X fini2 X fin3 and X finj4 in Eqs. All ID . (11191) . f II 2 3 j) and (I124p 
can be summed to give 
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In the last line we expanded around D = 4 and made use of Eq. ( fTBT) , 
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In Appendix C we evaluate the integral (11271) . The result (13121) . is a de Sitter invariant which 
does not contribute a leading logarithm. 
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We are now going to evaluate the one vertex integrals X sg n (n = 1, 2, .., 5). In order to do so, 
we shall need the basic integrals 
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(130) 
(131) 
(132) 
(133) 



The Y n and W n integrals are evaluated in Appendix B. In the asymptotic regime when a —>■ oo 
these integrals are given by Eqs. (12211 . (12221 . fl223|) and by fl224ti225l) . The S a integrals can be 
found in Appendix B3, Eqs. (f233ti243ll . while the A a integrals (a = D - 4, D - 3, D - 2, D - 1) 
are in Appendix B4, Eqs. (1215H21)UI) . 
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Recall that the single leg integral (15514551) is, 
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The singular part of this integral X^ sg can be combined with the singular one leg integrals X 



sg,l 



and X sgi3 in Eqs. (TTT2]) and (TTT9D to yield, 
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where (y+±) is given by Eq. (fl5|) . The terms in curly parentheses can be expanded in powers 
of D — 4. The result is, 



5 + (£>-4) -4 
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6 2 

This means that the most divergent contributions from the two left diagrams in Fig. 1 cancel. 
This is an example of a general phenomenon of cancellation of overlapping divergences and 
simplifies our calculation considerably. 

A further simplification is facilitated by the identity, 
dy dy 



-f(y) = -H% p (2-y)I[f](y) + V„V„/ 2 [/](y) 



(137) 



dx^ dx p ' 

where I[f](y) stands for the indefinite integral with respect to y. By making use of this identity 
and Eq. (11361) we can express the integral in f l 1 3 5 f) in terms of two relatively straightforward 
integrals as follows, 
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In order to accurately extract both the divergent and finite contributions from this integral, 
we shall keep the D dimensional form of the terms in K up to y 4 ~ D / 2 and y 2 ; in higher order 
terms y n (n > 3) we set D = 4, since these yield contributions that are suppressed linearly in 
D — A. Making use of Eq. ( jl"5l) we find, 
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from which it follows, 
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When these expressions ( I14UH1411) are inserted into Eq. f)138p and one makes use of the inte- 
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grals (I233H243I) evaluated in Appendix B3, one arrives at the following expression for X sgi A (11381) . 

V2 / v ; V 2 J (D-3)(D-4) 



55-81* 

(An) — 



x < 9n\p9t,]» 



2(D-1) 

2V(D-2) 



D-A 



1 \ „ [D 2 -1D- 



1 2(1) 2) D 1 J ~" 1 V 4(£>-3) D-2 



£>-4\ „ 



+ 



+ 



r(f)r(f+i)^ 

(D-2) 2 (Z>+8 
2(L>-4) 



D(D-2)(D-7) D 2 -7D- 



2{D + 2) 



12(L>-4) 2(D-3) 
r(D - 2) „ 



-D-4 



+ 



V[^][qVp] 

tf 2 



r(f)r(f +2 ) 

D-4 



(D-l)(D-2)"^ (D-2)(D-3) 



-D-3 



D 2 -7D + 8 „ 



+ 



8r(D - 2) 
(u-4)r(f)r(?+i) 



(142) 
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2(D-3)(D-4) 
Recalling that the operator, 

^2 ~ a %9u\[aO p] a — - g^ p g a]v a— 

Eq. (11421) can be evaluated and expanded in powers of D — 4. Here we provide an intermediate 
result, which also shows the terms linear in (D — 4), 
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The final result is, 
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The second part of the singular contribution to the one leg integral is the sum of Eqs. (11141) . 
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This can be expressed in terms of the integrals evaluated in Appendix B5 (as usually, the 
expression below is anti-symmetrized in both pairs of indices [[/,, v\ and [p, a]), 
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The result is, 
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Upon expanding this in powers of D — A one obtains, 
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What remains to calculate is the integral (I134p . This integral can be expressed in terms of 
the A a integrals ( 12441) . Firstly, we make use of Eqs. ( 11371) to write ( 11341) as 
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Secondly, we make use of Eqs. (11401 - 11411) and the definition of A a integrals (12441) to obtain (cf. 
Eqs. »), 
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This can be evaluated by making use of Eq. (11431) and Eqs. (I245H250"]) . Keeping the terms up to 
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order (D — 4) we have, 
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Upon rearranging the divergent contributions and dropping the 0(D — 4) terms one obtains, 
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Summing the integrals (JHQ, dM, ffl45|) . ffT49|) and (TT53D we finally get 
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This corresponds to the following values for the D-dependent quantities Ex, E 2 and F\ defined 
in expressions (j4"5]) and (I63ti64p . 
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n 

Our result for E 1 agrees with the stochastic prediction (l69il [2| , and our result for Fi agrees with 



the prediction (170|1 based upon conservation. 



VI. DISCUSSION 



By substituting (1155H157P into expressions (I65H66[) we obtain two loop results for the field 
strength contributions to the pressure p and to (p + p), 
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The M„) term in « agrees exactly with the stochastic prediction 2]. Together with the two 
loop results for the expectation values of (D^ipYDyip and (p*tp [l(, this constitutes an impressive 
verification of the stochastic formalism. 

Note that we have not bothered to work out the terms which fall off like powers of \ja. 
Just as for the stress-energy tensor of Ay? 4 — expressions (I50„5ip — these terms are separately 
conserved and play no role at late times. It has been conjectured that the 1/a contributions 
can be subsumed into a perturbative correction of the initial state 33[. Such a correction is in 
any case inevitable because free vacuum cannot be a very realistic initial state of an interacting 
theory. 

We have also not bothered to renormalize the composite operators F^F^, (D^(p)*D u (p and 
(p*(p. Recall that conventional renormalization only removes divergences from noncoincident 1PI 
functions. One would need to additionally subtract a series of counter-operators from F^F^ 



and the others 
finite part [62, 



;o remove their divergences, and there would be the usual ambiguities about the 



63). This is neither necessary nor desirable. The stochastic formalism which it 



was our purpose to check makes unique predictions for the leading logarithm contributions to 
the dimensionally regulated expectation values of the three composite operators [2J, and these 
predictions agree in each case. 

Accuracy is always an issue in such an intricate computation. Our result essentially consists of 
three numbers: E±, E2 and F±, which were defined in expressions ( 1451) and ( 1631I64T) and reported 
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in (I155H1571) . An obvious check on E% is that it agrees with the stochastic formalism for which a 
compelling and independent theoretical justification exists [2]. A powerful check on Fi is provided 
by partial conservation within the electromagnetic sector ( I6TH62I) . No direct check exists for E 2 
but it was of course computed using the same reduction strategy and many of the same integrals 
that produced correct results for E\ and F\. 

Combining the field strength contributions (11 5811159]) with the corresponding scalar results 
( 1751175]) gives the total for the stress-energy tensor of SQED (l49j) . 



Note that the divergent field strength and scalar contributions to (p + p) have canceled. This is 
required by the fact that the expectation value of T^ u gives the matter contributions to the 1- 
graviton 1PI function in gravity + SQED. Hence the expectation value of T Mi , can be renormalized 
by purely gravitational counterterms, and all of these degenerate to constants times g^ on de 



see that expression (11601) for p can contain divergent constants, but the ln(a) terms in p and all 
of {p + p) must be finite. Because the constant contribution to the pressure can be absorbed into 
a renormalization of the cosmological constant, this term has no physical significance — which 
is reassuring because it is the least well checked. 

Although this exercise was undertaken to check the stochastic formalism, our results (11601 - 
11611) have considerable physical interest in their own right. After many e-foldings the In (a) 
contribution to the pressure fll60j) must dominate and we see that SQED induces a growing, 
negative vacuum energy. The physical interpretation is that the vacuum becomes polarized by 
the inflationary production of charged scalars. Just as a dielectric slab will be drawn into the 
region between the plates of a charged capacitor, so the production of additional scalars seems 
to be favored by the electric fields of those that came before. 

It is natural to wonder how far this progresses, both at higher orders in the loop expansion 
and in time. Recall that the two loop ln(a) derived entirely from the field strength contributions. 
That is not true at higher orders, however, only the field strength contributions act to decrease 




Sitter background. Because the expectation value of T^ u takes the form pg^ u + (p + p)a 




33 



the vacuum energy [$|. At £ loop order one gets £ — 1 factors of e 2 /(47r) 2 , and there can be 
up to £ — 1 factors of ln(a) The leading logarithms at all loops become order one after 
ln(a) ~ (4-7r) 2 /e 2 e-foldings. At this point perturbation theory breaks down and one must employ 
a nonperturbative resummation scheme to evolve further. It was to solve this sort of problem 
that the stochastic formulation of SQED was developed, and the answer is known [2J], 

A 

lim pit) « 0.6551 x — — « — — x 0.2085 x GH 2 . (162) 

Here A = 3H 2 is the bare cosmological constant and G is Newton's constant. 

Although (11621) is a nonperturbatively large decrease in the vacuum energy, it is suppressed 
by GH 2 relative to the vacuum energy of the bare cosmological constant. The largest value 
of GH 2 consistent with the current upper limit on the tensor-to-scalar ratio for anisotropies in 
the cosmic microwave background is about 10 12 [66], so ( 11621) does not represent a significant 
decrease of the vacuum energy. On the other hand, it is an enormous amount of vacuum energy 
by current scales. Further, the shift is dynamic; it was caused by inflationary particle production 
and it would presumably dissipate, on some time scale, after the end of inflation. This may well 
have important consequences for cosmology [23 ] . 
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Appendix A: Extracting d'Alembertians 

Here we list the results of extracting the d'Alembertians from powers of y/4 which reduce the 
degree of divergence by 2. 

When a d'Alembertian acts on a nonsingular function F = F(y) one obtains (11061) . 

^-F(y) = (4 - y)yF"(y) + D(2 - y)F\y) . (163) 
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(A nonsingular function F of y — y(x; x') is a function which, when expanded in powers of y 
does not contain the power y 1 ^ / 2 .) Equation (11631) is helpful for establishing the relations 
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which when applied to /3 = (D/2) - 1, (D/2) - 2, (D/2) - 3, (D/2) - 4, -2, -1 yields 
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Appendix B: Some useful integrals 



Bl. The Y integrals 



Here we evaluate the basic integrals of the form, 
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In order to perform the integrations, the following relations are needed, 
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We shall also need the polylogarithm functions, defined as, 

U n (z) = PolyLog[n, z\ = — (n = 2, 3) . 
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Note that the natural logarithm is a special case of the Li function, ln(l— z) = — Lii(z). Moreover, 
we shall make use of the properties of polylogarithms, 
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The simplest nontrivial integral (n — 1) in Eq. (11931) is, 
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where t)q denotes an initial conformal time for which we choose a(r] ) = 1 and we made use of 
the surface area of the D — 2 dimensional sphere, 

S^ = ^L = ^p^l. (205) 
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The integral (12041) can be expressed in terms of a hypergeometric function. Since we are ulti- 
mately interested in expanding in powers of D — 4, it is more convenient to expand the integrand 
in powers of D — 4, 
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and then to integrate. The result is, 
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The next integral in Eq. (11931) is y 2 - With the help of Eqs. (TT951) . (1206]) and (12751) we can 
write (x = r/ Ai], v = a' /a, u = 1 — v), 
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where we introduced a new variable, u — 1 — a' /a. These integrals can be integrated to give, 
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The last integral in Eq. (11931) is Y%. Similarly as in the case of Y2, we take account of Eq. f l 1 9 6 f) 
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Upon performing the integrals in (12 10[) we finally obtain 
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B2. The VP integrals 



The following class of basic integrals we need (these integrals we calculate only up to order 
(D — 4) 2 ) is of the form, 
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The trivial integral is Wq = 0. The first nontrivial integral can be written as 
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After the radial integration is performed, this integral can be reduced to, 
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The second integral Wi in (12121) we only need to leading order in the D — 4 expansion. We 
begin with, 
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After the radial integration this integral becomes, 
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We shall now act with d'Alembertians on the basic integrals (j207p . (I209p and (12111) . For 
simplicity we neglect the terms that are suppressed as In 2 (a) /a, or more since they vanish in the 
limit when a — > 00. Since Y n are functions of the scale factor only, the d'Alembertian operator 
simplifies to, 
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we can use the asymptotic form of the basic integrals (I204p . (12091) and (I21ip ) to calculate the 
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late time behavior. In the limit when a — > oo, the integrals become, 
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The asymptotic forms for the corresponding (nonvanishing) d'Alembertians are 
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and similarly for W n '8. 



B3. The 5 integrals 

We are now ready to calculate intermediate integrals that we need for evaluation of the one 
vertex integrals. They are of the general form, 
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with a = ^ - 6, ^ - 5, ^ - 4, 3 - ^, D - 2, D - 3, D - 4, § - 1 and § - 2. All of 5 Q integrals 
can be represented in terms of the basic integrals Y n given in Eqs. ( 122 11 - 12231) and the associated 
d'Alembertians (1226H23TD . 
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D 



-2, 



'—2 



a x a 



■D-4 
2 



ip - 4)y x + ±p - 4) 2 F 2 - ±p - 4) 3 F 3 + O (p - 4) 4 ) 



rp) 

+ P-4) 3 



p-4) 



1, , , 11 

— In (a) H 

2 v ; 12 



P-4) 5 



91 



-In 2 (a) ln(a) H 

8 w 3 w 144 



_ J_ i n 3( a ) + I ln 2( a ) + ( _I _ H ] ln ( a ) 
48 V ; 6 V ; V 6 48 ' 



+ 



257 llvr' 



1728 288 



+ 0(p-4) 4 ) j. 



(234) 



44 



The following integral is H3d__ 5 , which we can rewrite with the help of Eq. fl 168[) as, 
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Next integral we calculate is a = y — 1. With the help of Eq. (11741) . we arrive at, 
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j In (a 



/ 89 


7tt 2 V 


V 192 ^ 


96 ) 



+ 0((D-A) 4 ) j 
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In order to evaluate the integral a = ^ — 4, one needs to subtract the divergent contribution 



according to Eq. (11671) . The result is, 

2 □ D-8 



+ 



(3D - 10)(D - A) H 2 2(D-A)_ 
□ D(D-2) 



(237) 



(3D-10)(D-4) [ H 2 
,(4vr)fr(f) 



4 



(4tt 



iJ 15 (3D-10)(D-4)r(f-l) 



# D T(D) 



(D-Af 



13 



12 + p-4) 



65 



4 + 



7T 



21 



— In 6 (a) ln z (a )+[ 



61 7tt' 



6 



ln(a) 



19 13tt^ 



4 



+ C(3) 



+ 



T(D) 



r(f)r(f-i) (3D-io)p-4) 



+ O ((£> - 4) 3 ) I . 



(238) 
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Furthermore, for a = 30 - 3 W e make use of Eq. (I166p to obtain, 



□ 



D - 6 



(3D - 8)p - 3) H 2 2p-3) 

Wj£(f) J 



# D rp)\ 

+ p-4) 2 



12 + (D-4) 



14-1- 
2 



_ l^l n 3(a) - — ln 2 (a) + 
4 v ; 2 v ; 



107 7tt 2 

"2T + ~Y 



ln(a) 



rp) 



115 5tt 2 
— + - - 8C(3) 

D- 6 



+ 0(P-4)»)j. 



r(f)r(f-i) p-3)(3D-io)p-4) 

The next class of integrals begins with a = D — 4. The integral is of the form, 



(239) 



— D-4 



-p - 4)n + ~p - 4) 2 F 2 - ~p - 4) 3 F 3 + o (P - 4) 4 ) 
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+ p-4) 3 

+ 



(13-4) 
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In 2 (a) - -ln(a) + 



35 _ tt 2 
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-ln 3 (a) + -ln 2 (a) 



32 tt 2 " , 
T +12) Ma) 



/217 


7tt 2 V 
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24 J 



0{(D-4f) 
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With the help of Eq. (11721) we calculate the next integral, 
2 □ 6 



— D-3 



= I 



(D-4)(D-6)H 2 D-6 
; (4vr)fr(f) 



-D-4 



H D T(D) 
(D-4) 3 



3-4p-4) + (D-4)' 
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iln 3 (a) + ^ln 2 (a) 



7T 



x ) ln (°) 



llTT' 

~Y1 



(241) 

0(p-4) 4 )}. 



Next integral corresponds to a = D = 2, and can be calculated by the means of Eq. (11711) . The 
result is, 

2 □ 4 



— D-2 



+ 



+ 



(D-3)(D-A)H 2 D-4 
2 [ □ D(D-2) 

(D-3) (D-4) [IP + 4 

, W#r(f) 
# D rp) 

+ p-4) 2 

rp) 



+ 



(47T) 



p-3)p-4)r(f-i) # D 



12+ p-4) 



4 + 



7T 



Il n 3( a )_^ ln 2 (a) + 

4 V ; 8 K J 



1 
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5vr 2 \ , . , /25 53vr 2 
ln(o) + — + 



4 
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24 



r(f)r(f-i) p-3)p-4) 



+ G(p-4) 3 ) I 



(242) 
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The final integral corresponds to a = D — 1 and it evaluates to, 

2D 2 



— D-2 



(D-2) 2 H 2 D-2 

(47r)fr(f)r 



7T 



+ (D-4) 2 

T(D) 



7, o. , 7, 2 . . /193 5tt 2 \ , . , /37 vr 2 

- ln 3 (a) n 2 (a) + 1 ln(a) + 

4 v ; 2 v ; V 24 4 / v ; 1 24 12 



T(f)r(f-1) (D-2)(D-3)(D-A) 



0((D-A) 3 ) 
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B4. The A integrals 



There is one more class of integrals to be evaluated which are based on the basic integrals (12X21) 
that are needed to evaluate the one vertex integral contributing to the central diagram in Fig. 1. 
The general form of these integrals is, 



A, 



d D x'a ,D ln(a') 



2/4 



4 



2/H 



4 



(244) 



with a = D — 1, D — 2, D — 3, D — 4, y — 1 and -f — 2. All of these integrals can be represented 
in terms of the W n functions given in Eqs. (12151) and (12181) . We give these integrals to order 
(D — 4) 2 , which is also what we need for the evaluation of the integral (11 34ft . 

When expanded in powers of D — 4, the K a=D _ A integral can be written as, yields 



A 



D-4 



d D x'a' D ln(a') 



2/+H 



-0D-4) 



4 



-{D - A)W 1 + - {D - A) 2 W 2 + 0({D- 4) 3 ) 



, W#r(f) 
h d r (d) 



(D-4) 



+ {D-Af 



4 



l l2n 11 

-In 2 (a) H 

2 v ; 6 



35 



ln(a) 



85 
36 



(245) 
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7T 



- ln 3 (a) - - ln 2 (a) + — - — ln(a) - — + 



9 



6 



125 llvr 2 



27 



36 



C(3) 



+0({D-4) 3 ) 



Similarly, for other integrals we need we have, 



A, 



D -4, 



-Wi 



(D-±y 



W 2 + ((D - 4) 3 ) 



.(4 7 r)fr(f) 



H D T(D) 
' 1 
24 



(D-4) 



--ln 2 (a) + — ln(a) - 
4 v ; 12 v ; 72 



+ {D-4f 



In 3 (a) 



In 2 (a) 
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(246) 
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+ 0({D-4) 3 ) 
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A 



D-3 



□ 



6 



(D-4)(D-6) H 2 (D — Q) 



i7 D r(D)\ 



+ (£>-4) ; 



-31n(a) + 
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9 



+ (£>-4) 



-41n(a) + 6- 



7T 



7T 



- ln 3 (a) - - ln 2 (a) + 7 - — ln(a) - 7 + — + 3C(3) 



7T 



(247) 
+ 0((D-A) 3 ) 



A £ _ 1 
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^{-3 ln(a > + >- 4 > 

+ (/J-4) 2 \\\n\a) -\\n\a) 
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49 
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.(4tt)t 1 2 
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; (4vr)fr(f) 
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r(D) 



_ll n 3( a ) + ^i ln 2 (a) 

24 V ; 12 V ; Ml 



2323, , . 2533 9 
ln(a) + — - 2vr 2 
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r(f)r(f-i) (/j-3)(/j 
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(249) 



Ad_i = 



□ 



A^ 2 + C? {(D — 4) 3 ) 
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(D-2) 2 H 2 D-2 



+ P-4) 
+ P-4) 2 

r(D) 



| ln 3 (a) + 3 ln 2 (a) + (-g - 2vr 2 ) ln(a) + ^ - 3vr 2 + 12C(3) 

.1! i n 3 (a) _ 361 ln2(a) / 169 2 \ 233 2 _ 
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r(f) 2 ( J D-3)(/j-4) w r(f) 2 P-2)P-3)(/J-4) 



+ 0((D-A) 3 ) 



(250) 



48 



B5. The 0, integrals 



We can now use these results to calculate the following class of integrals we need, which are 
of the form, 



<x,P 



a x a 



a/3 



o. 



4 



-a. 



0, 



?7h 



(a + p)(a + /3+l) 



V„V 



afiH 2 g iip 



P P 



"v 4 
1 



-0 



(251) 



□ 



D - 2(a + /3 + 1) V (a + /?)(a + /? + 1) # 2 



where we made use of Eq. (I184p and assumed a+/3+l ^ . These integrals can be evaluated with 
a help of Eqs. (I185H192]) and the H 7 integrals (I233H242]) . We need the integrals corresponding 
to the pairs, (a, p) = (f - 1, D - 2), (f - 1, D - 2), (f - 1, D - 3), (f - 1, D - 4), (f - 1, D - 
5),(f -l,D-6),(f -2,D-2),(f -2,D-3),(f -2,D-4),(f -5,D-5),(f - 1, f - 



2),(f -3),(f -2,f -2). 

Note that when V M V P acts on a function of time (or equivalently of a, which is a scalar 
function), one obtains, 



V^Vp 
H 2 
□ 

H 2 



-9, p ad a + (a 2 ft°) a 2 d 2 
gvYjjV* = -a 2 ~ D d a a D d a . 



(252) 
(253) 



Making use of these relations and of Eq. (11841) (and assuming that a + (3 + 1 ^ ~), Eq. (125 ip 
reduces to, 

a 2 ~ D d a a D d a 



n 



H z g, P 



+ H 2 (a 2 6y p ) 



a(3 

a+p){a+P+l) 

a(3 



ad a + 



a(3 



D-2(a+f3+l) \(a+(3)(a+(3+l) 



(254) 



(a + /3)(a + /3+l) 

Note that the last term breaks the Lorentz symmetry. The £l a ,p integrals are symmetric under 
a <=± p. 

Since the one vertex integrals are at least linearly divergent in D — 4, to extract the finite 
contribution, we only need terms up to linear order in the D — 4 expansion. 

We begin with the integrals of the form a = ^ — n, P = D — n (a + P — W- — (m + n); m = 
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1,2, 3, 4; n = 2,3,4, 5). We get, 
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The following integrals are also useful, 
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3 v ' 9 3 

Next class of the f2 integrals is of the form, a = -j — m, (3 = y — n (a+f3 = D — (n + m); m,n = 
1, 2, 3). Since the coefficient of these integrals is not singular when D — 4, actually we need only 
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the 0((D — 4)°) contribution. The integrals are, 
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B6. The radial integrals 



The radial integrals are of the form 



J n = f dxx D ~ 2 [\n(l - x 2 )] n 
Jo 



and when evaluated result in 
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Appendix C: The finite two leg integral 



In this appendix we evaluate the two leg finite integral (1127p . The relevant part of the integral 
is of the form, 



J MP = ^(±)(±) / dVa' 4 f d A x"t 



7/4 



d 1 







In 



/V±± 



(278) 



v cta^ y + ± J \ dxP y'\_± J \ 4 
The integral can be naturally split into the spatial and timelike parts, each of which we evaluate 
separately. The spatial part of (12781) can be written as, 



T. 



3a 2 H 4 



d 4 x'a' 3 f d A x"a"* . ja'a"H , 2 
Af • Af " In ( : Ax' ±± 
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(279) 



It is natural to introduce the coordinates r 



\x — x'\\, r" 



\x — x 



in terms of which 



\x' — £"11" = r 2 + r" — 2rr" cos(#) , where # = Z(x — x', x — x"). The angular integral over 9 
can be easily performed. The result is, 
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4 
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/A 2 

2 



Ar/" 
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^280) 



where we also used the symmetry of the resulting integrand under the exchange r — > — r to 
extend the limit of integration of the r integral from [0, oo) to (— oo, oo). Finally we made use 
of the fact that temporal integrals extend from t]q = — 1/H to r/, and that the symmetry allows 
us to constrain the rj' integral up to r/', thus gaining a factor 2. With this we have achieved 
an important time ordering, rj' < rj < rj, which we will find useful below. Next we split the 
integral (I280p into two integrals lij = 2^- + 2^- as follows, 
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(282) 

(283) 
(284) 
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X 


X 


X 


Ar/ 


A 

IA7] 


A 
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+ 


+ 


+ 


— ie 


—ie 


—ie 


+ 


+ 




—ie 


+ie 


+ie 


+ 




+ 


+ie 


—ie 


—ie 


+ 






+ie 


+ie 


+ie 



TABLE I: The ie prescriptions of the intervals Ar/ = rj — rj , Ar/' = 77 — 77" and At/ = rj' — rj" as a 
function of the signature ± of the vertices x, x' and x". 



Before we begin evaluation of these integrals, we recall the ie prescription 

^ 2 » ||Ax|| 2 - (A V -ie) 2 



Axl 
Ax 2 + _ 
Ax 2 



Ax\\ 2 - (\Arj\ - ie) 
Ax\\ 2 - (Ar] + ie) 2 , Ax 2 _ + = || Ax|| 2 - (Ar/ - ie) 
Ax\\ 2 - (\Ar)\ + ie) 2 -> || Ax|| 2 - (Ar/ + re) 2 , 



(285) 



where the implication follows as a result of the time ordering 77 > rj' > rj" in Eqs. (I281H2821) . 
Analogous prescriptions hold for the other two intervals. Based on this, we arrive at the ie 
prescriptions as given in Table I. Note that Ar/' and At/ have identical ie prescriptions and 
that they are completely specified by the signature of the x" vertex. On the other hand, the ie 
prescriptions of Ar/ are completely given in terms of the signature of the x' vertex. 

We shall now use this fact to perform the r integrals (I283H2841) . The first integral JW can be 
performed by making use of the Dirac identity, 

r r 1 
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d d 
—5(r + Ar]) + — 5(r - Ar/) 
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[r 2 - (Ar/ - ie) 2 ] 2 [r 2 - (Ar/ + ie) 2 } 2 4Ar/ 

The derivative d/dr can be moved by partial integration to the integrand, such that the result 
of integration of (12831) is, 
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. (287) 
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Analogously, to evaluate we can use the following Dirac identity, 
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with whose help one can evaluate J( 2 \ The result is, 
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The symmetry of the integrals and J*- 2 ** under r" — > —r" allows one to extend the r" 
integration to — oo. Upon inserting these two integrals into Eqs. (1281M2821) and summing them, 
we obtain, 
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It is convenient to rewrite this expression in terms of two radial integrals as follows, 
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can be partially integrated to give, 
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The first term in can be rewritten as 
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resulting in four simple integrals, + = La + Lb + Lc + where 
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where we made use of r" + Ar] = r" — Ar/" + Ar/ + Ar/" = r" + Ar/" — Ar/'. All of these integrals 
are simple to evaluate by contour integration. For example, upon writing the ie prescription 
explicitly (e.g. from Table I), La can be written as 
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This is straightforward to evaluate, 
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Similarly, the remaining integrals in Eqs. (I297H299"]) evaluate to, 
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Upon combining these integrals, we can write the spatial integral Z^- (12911) as follows, 
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Upon changing the variables to v = a' /a and w = a" /a, this reduces to 



Z 



ij 



32tt 4 a 2 5 tJ [* r 

' aw / dww 



3 H & 



l/a 



l/a 



In 1 



— ln(l — w) — 
vl v ; 2v(l-w) 



(306) 



These integrals are straightforward to do. The result is, 
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This is the final result for the spatial part of the finite two vertex integral. 

The de Sitter invariant form of the integrand in Eq. (I278j) indicates that the leading order 
contribution to the integral X^ u should be de Sitter invariant, such that when combined with 

4 

the 00 component of one expects to get, = j^g^ plus terms suppressed by powers of 
the scale factor a. A detailed evaluation confirms this expectation. Albeit more cumbersome, 
the evaluation of Zoo closely resembles the evaluation of Zy, and here we present only the major 
steps. 



From Eq. (12781) it follows 
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Upon performing the r integral one obtains, 
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Performing the r" integral requires more work, 
32vr 4 rv 



Z 



00 



tf 2 



no 



drj'a' 3 / drj"a" 3 



Ar/Ar/'ln a' a" H 2 (Ar]") 2 - Ar]{Ar)" + Ar/) 



+ 77 



+ 77 



(A77 + Ar/') In a'a"H 2 (Ar]"" 2 



2Ar}' 



In 



a'a"# 2 Ar/'Ar/ 



1 A77 



2 Ar/" 



(310) 



57 



By making the substitutions v = a'/a and w = a" /a, this can be reduced to a set of relatively 
simple integrals. The result of integration is, 




When the two parts 2^- in Eq. (I307j) and X Q0 in Eq. (13 lip are combined, the leading order 
contribution is de Sitter invariant, 

J ^ = Jje^u (312) 

plus terms suppressed as powers of scale factors, in accordance with the expectation. 

From this result we see that the finite integral Zfj n does not contribute a leading logarithm to 
the two loop photon field strength bilinear. 
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